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ABSTRACT 

The generation and the propagation of three fundamental tidal modes - 
the symmetric diurnal mode , the symmetric semidiurnal mode and the 
antisymmetric diurnal mode - have been studied based on a theory outlined in 
a foregoing paper. The equivalent depths of these modes and the eigenfunctions 
which describe the latitudinal structure in the pressure field are shown to change 
wuth altitude. While within the lower atmosphere the eigenfunctions are the Hough- 
function they transfer into the ordinary spherical functions at thermospheric 
heights. The transition occurs between 100 and 200 km. Horizontal winds 
and pressure amplitudes of the three modes which depend on altitude and lati- 
tude are determined and compared with available wind and density data. 


INTRODUCTION 


A. 

In the first part of this paper (Volland and Mayr, 1971a; referred to as 
paper I) we developed a theory of thermospheric dynamics in terms of the 
eigenfunctions of the dissipative atmosphere. We found approximate general 
analytic solutions for the generation and propagation of tidal and planetary 
waves excited by the solar heat input due to XUV-radiation and corpuscular 
heating. That solar heat input has also been developed in terms of eigen- 
functions and contains tidal components (depending on local time) and planetary 
components (depending on seasonal time). For a quantitative evaluation of 
the propagation characistics of the various atmospheric wave modes the eigen- 
values of the waves must be known. These eigenvalues depend on the equivalent 
depths which in turn are responsible for the latitudinal structure of the corre- 
sponding eigenfunctions. It is the purpose of this paper to calculate the equiv- 
alent depths of some of the important tidal modes on an approximate analytic 
basis and to discuss their height dependence. Planetary waves will be considered 
in the third part of this paper. A numberical calculation taking account of a more 
sophisticated thermospheric model will be presented in an additional paper 
(Volland and Mayr, 1971b; referred to as paper II). 

2. APPROXIMATE ANALYTIC EVALUATION OF THE EQUIVALENT DEPTHS 

In section 4 of this paper I we determined an analytic solution (Equ. (1/32) 
for the wave amplitudes of the mode with eigenfunction 
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Here, P™ are the associated Legendre polynomials in Schmidt's 

normalization, (r, A ) are the spherical co-ordinates, z=r-R o is the 

height above the ground, f = «/2 n is the Coriolis parameter, cc is the 

angular frequency of the wave mode, 12 is the frequency of the sidereal day, 

12 is the frequency of the sidereal year and t is the universal time. 

The solution (l/32) can be evaluated provided the eigenvalue and the 

corresponding equivalent depth of the wave mode are known. The equivalent 
in f 

depth h R * was defined in the usual way (Chapman and Lindzen, 1970) as the 
solution of the eigenvalue problem of the following set of equations: (see 
Equs. (1/16) and (1/19) ): 
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are the relative wave amplitudes of pressure, northerly wind and westerly 
wind, and 
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are series of sperical harmonics to be determined below. p Q is the mean 
pressure of the steady state thermosphere, c q is the velocity of soundy~1.5 the 
ratio between the specific heats , 
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Z yis an d Z cq1 are the disspiation factors of molecular viscosity and ion-neutral 
drag (see Equs. (1/3) and (1/4), and g is the earth's acceleration constant. 

The latitude dependent part of the eigenfunction in (1),0™’^, degenerates 
to the well known Hough -function within the nondissipative atmosphere (Chap- 
man and Lindzen, 1970). It degenerates to the spherical surface if the 
Coriolis force can be neglected (see Equ. (1/23). That is, 


surface function P™ if the Coriolis force can be neglected (see Equ. (1/23) ). 
That is, 
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Because of Equ. (6), we shall adopt the following definition for the 

in f 

numbering and the normalization of the eigenfunction 0™* : the wave 
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domain numbers (n, m) of the eigenfunction 0 shall pass over into the 


wave domain numbers (n, m) of the corresponding spherical function P. 
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Then, in the case of a highly dissipative atmosphere, the function 0 


n 


approaches the form (6). It should be mentioned that our numbering and 
normalization differs from other definitions of the tidal Hough -function 
(e. g. , Chapman and Lindzen, 1970; Tarpley, 1970). 

TX1 f 

An exact numerical determination of the equivalent depth h™’ 1 can be 
found from Equs. (3) applying the usual expansion procedure (Chapman and 
Lindzen, 1970) though our system is mathematically somewhat more diffucult 
due to the dissipation factots Z. Within the thermosphere the dissipation 
factors Z increase with altitude and ultimately dominate the Coriolis force 
if | Z| » 1. That is generally the case above 200 Inn altitude where the eigen- 
function 0^ ’ in fact approaches rather well the associated Legendre poly- 
nomial pJJ 1 (see Fig. 1). Therefore, we shall adopt in this paper an approximate 
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method to find from the systern(3) analytic solutions of the equivalent depths 
which sufficiently well converge to the corresponding values for the non- 
dissipative atmosphere or to the formulae of the equivalent depths of Equ. 
(1/26) in the case of a highly dissipative atmosphere. We do this in the 
following manner: We apply the assumption for the eigenfunctions of Equ. 
(1) with the normalization from Equ. (6). We retain the two first equations 
of (3), but replace the third equation of (3) by 
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Since the Hough-functions are tabulated, we can compare those exact 
data with our approximate calculations in the case of zero dissipation. It will 
turn out that the Hough -functions are sufficiently well represented in our 
procedure while the equivalent depths differ by not more than 20% when 
compared with the exact values. Moreover, our wind fields, though very 
similar in the latitudinal structure, differ somewhat from the exact structure. 

In view of the advantages of simple analytic solutions which allow to discuss 
carefully the dependence of the eigenfunctions on heights, frequency and 
dissipation we shall, however, tolerate these differences. The exact treat- 
ment will certainly modify our approximated results, however, it will not 
change their substance. We should add that this approach is valid only for 
the lowest zonal wave domain numbers n. 

In the next section we shall calculate explicitly the eigenfunctions and the 
equivalent depths for some of the important tidal modes. There we shall restrict 
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ourselves to those tidal inodes which appear to be generated most 
intensively at thermospheric heights (see Eq. (1/9) ) and which form the 
predominant terms in the series of spherical functions for the observed 


temperature distribution (see Eq. (1/14) ). These tidal waves are (in our 


notation) the symmetric diurnal mode 0 


1 , 1/2 


, the antisymmetric diurnal mode 


_1 j /2 2 1 

y 2 * and the symmetric semidiurnal mode 0 2 * which degenerate at 

11 2 

thermospheric heights into the spherical functions P^ , P 2 and P 2 . We 
shall study the change of their latitudinal structure with height and their 
propagation characteristics. 


3. SYMMETRIC DIURNAL TIDAL MODE 0 

3a Eigenfunction and quivalent depth 

The method to solve the eigenvalue problem of Equ. (3) has been described 

in some detail for the fundamential diurnal mode called the (1, -1) mode by 

Kato (1966) or the (1, -2) mode by Lindzen (Chapman and Lindzen, 1970). 

That mode has been discovered independently by Kato (1966) and Lindzen 

(1966). It is a trapped mode within the lower atmosphere and the main 

generator for the geomagnetic Sq current (Stening, 1969; Tarpley, 1970). It 

plays a dominant role within the thermosphere (Volland and Mayr, 1970). 

1 l/2 

We shall see in the following that in our notation this mode is the 0 ^ 

-mode. 

From the tables of the Hough-functions (Chapman and Liudzen, 1970) we 
find an expension into spherical functions of the diurnal (1, -2) mode of 

0(1, -2) = 0.777 Pj + 0.582 P* + 0.0751 P* + . . . 
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where the neglected terms are of the order of 5% and smaller. In order 
to approximate this Hough -function we start from an expension. 

-Mi-; 


(9) 


where we already used our numbering and normalization cheme from Eq. (6). 
For f = 1/2 (solar diurnal component) and Z = 0 (nondissipative atmosphere), our 
eigenfunction (9) should converge sufficiently well to Eq. (8). We futher- 
more assume for the series of the horizontal winds in (4) 
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in (3). This restriction will be abandoned in the numerical treatment of 
paper II. 

Substituting Eqs. (9) and (10) into the first two Eqs. (3), we find after 


some calculation the exact solution 
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Here the numbers bn the r.h.s. of (12) have been determined for f^=f = l/2 
Introducing the velocity components of Equ. (12) into Equ. (7) we notice 
that the r.h. s. of (7) is fulfilled and that the equivalent depth is 


l,f 2fl 2 r 2 f A 
1 = g(f k -l/6) 


(13) 


In order to compare our result with the corresponding Hough-function 

0 (1-?.) of (8), we set f ^ = f = 1/2. Then we obtain from (12) 

*1,1/2 _ 8 V5 - _ 

°1,3 --27 - - 0 ' 725 

which must be compared with the corresponding number from ( 8) of 


0.582 

0.777 


0.749 . 


Both numbers differ by not more than 3%. The latitudinal structure of 
both representations of the Hough -function (1, -2) are plotted on the bottom of 
Fig. la as solid and as dash-dotted lines indicating the excellent approximation 
of our function. 

The exact numerical value of the equivalent depth of the(l,-2) mode is 
h(l, -2) = - 12.27 km (Chapman and Lindzen, 1970). 

We compare that number with our value from (13) 

h? - ’ ~ — = - 9. 9 km for r = R and f . = f = 1/2 

1 20 g ok 


which is about 20% smaller than the exact value. That difference could be 
compensated in principle by a virtual earth’s radius R q > R q . 
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The bottom curves in Figs, lb and lc show the horizontal wind fields 
,1/2 ^ ( j ) 1,1/2 our ap p rox i ma tion in the case = f = 1/2 (solid lines) 

and of the exact treatment (dash-dotted lines) thus indicating rather similar 
latitudinal structures. However, the differences are more pronounced than 
in the case of the Hough-functions of Fig. la. 

Within the dissipative thermosphere it is If, I » 1. Thus, 


6 1,f « 
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00 
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and the wave amplitudes of pressure and horizontal winds indeed approach the 
formulae of (1/20) to (1/26) where the eigenfunctions degenerate to the ordinary 
spherical functions. 


3b Wind and pressure fields 

In order to study the behavior of the eigenfunction within a dissipative 
atmosphere, we considered a model of the steady state thermosphere of 
Jacchia (1964) with an exospheric temperature of T ^ = 1000° K, and we 
assumed plausible data of the ion-neutral collision number from a paper of 
Dalgarno (1964) and dissipation factors Z v . g and Z^ from full wave calcula- 
tions of thermospheric dynamics by Volland and Mayr (1970). Th* tem- 
perature profile and the dissipation factors used are plotted in Fig. 2 versus 
altitude. Here it must be mentioned that due to our approximations in con- 
nection with molecular viscosity, our dissipation factor Z^ = % co i + z vig 
becomes ill defined above about 300 km altitude. 
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Below that height, it is generally Z qo j > Z ^ . However, the term Z j 

decreases with height above 300 km due to the decrease of the electron 

density. Between 300 and 400 km, the term Z cq ^ may therefore fall below 

the viscosity term Z . . On the other hand, due to this decrease of Z, . , 

the horizontal velocity may increase which in turn should lead to an increase 

of the viscosity term Z . . For convenience, we therefore assume that Z, . 

vis kin 

remains constant above 300 km altitude. Moreover, we consider an isothermal 
atmosphere below 100 km. 

1 i/o 

In Fig. la we plotted the eigenfunction© ' versus co-latitude d for 

the four different heights z = 100 km, 150 km, 200 km and 300 km in amplitude 

(solid lines) and phase (dashed lines). The phase is given in terms of the delay 

1 1/2 

with respect to p^* (in hours) according to 
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We observe in Fig. la a transition from the Hough -function at 100 km to 
the spherical function P* at 300 km altitude. That transition occurs mainly 
between 100 and 200 km. The sperical function P* is an excellent approxima- 


> 1 , 1/2 
1 


tion of 0*’ " at 300 km with an accuracy of better than 2%. Similary, in 

1 1/2 

Figs, lc we note a transition in the longitudinal wind function from a 

wind changing sign at d = 65° co-latitude at 100 km to a nearly constant 
wind velocity at 300 km which agrees with the assymptotic value of 




1 , 1/2 


® Pj/ sind = 1 
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in Eq. (1/10). The latitudinal wind function in Fig. lb remains 

unchanged with height in our approach and is proportional to 


dpl/di?= costf 


However, the amplitudes and phases (here in terms of the time delay 

1 1/2 

with respect to the pressure p^’ ) change drastically with height due to the 

increasing influence of the dissipation factor . 

1 1/2 

The equivalent depth h^ ’ of Equ. (13) is a measure of the propagation 
characteristics of the wave mode. Since the equivalent depth has no direct 
physical meaning we better describe this behavior by the eigenvalue q from 
Eq. (1/37), the real and the imaginary term of which are responsible for 
vertical phase velocity and attenuation. The height dependence of an up going 
wave is given by (see Eq. (1/33)) 


exp -{(j^ + 0 - 1) z/2H q J. . 

The vertical phase velocity can be represented by a vertical wave 
length of 


X 
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47T 

a 
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In Fig. 3 we plotted the attenuation factor ($ - 1) and the propagation 
1 12 

factor O' of the 0 * - mode versus altitude. The upper scale in Fig. 3b 

gives the normalized vertical wavelength X /H according to Eq. (15). From 

z o 

ll/2 

Fig. 3 we notice that the 0^’ - mode, while purely evanescent below 

100 km (Oi = 0), changes into a quasi - evanescent mode within the thermosphere. 
Here, wave attenuation reaches a maximum near 140 km and then diminishes 
below the 100 km -value within the upper part of the thermosphere. 
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The vertical wavelength is infinite below 100 km. It remains very large 
(of the order of 1000 km) at thermospheric heights. 

From Figs. 1 to 3 we conclude that the thermosphere above about 200 
km height as well as the lower atmosphere below 100 km behave nearly like 
isothermal atmospheres with constant parameters. Therefore our solutions 
derived in (1/41) and (1/50) can be applied to a sufficient degree of approxima- 
tion, while the transition region between 100 and 200 km must be treated by 
full wave theory to be discussed in paper II. 

1 1/2 

We calculated from (1/48) the generation and propagation of the 0 ’ 

-mode within an isothermal nondissipative lower atmosphere assuming a 
heat input between the ground and 100 km proportional to the mean pressure: 

1,1/2 (-0.01 for L z 100 km 

= (16) 

1 Po (0 for > z 100 km 

with the maximum heat input during local noon. We plotted the relative ampli- 
tudes of vertical wind w/c q , pressure p/p q and temperature T/T q versus alti- 

tude in Figure 4. Though theO^’ -mode is evanescent #3-1 = 0.25) in that 
height range, the integrated sum of the up- and downgoing waves coherently 
generated at any height interval below 100 ’nn add together to a net amplitude 
asymptotically approaching constant values for vertical wind and pressure 
(see Eq. (l/52b). On the other hand, the temperature decreases, except in the 
vicinity of 100 km which is the upper boundary of the heat input. 
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Above 100 km the wave amplitudes decrease like free internal wave 
according to Eq. (I/50c). This indicates the small penetration depth of 
evanescent waves which are* only significant within the altitude range of 

2H 

Az = ^ ~ 50 km ^ ^ 

in the case of upgoing free internal waves within the nondissipative lower 
atmosphere. 

The discontinuity of the temperature amplitude at 100 km in Fig. 4a is 
the result of the abrupt termination of the heat input at this height level in 
our model. Below 100 km, the temperature amplitude is generated by the heat 
input and by adiabatic expansion and heat advection due to winds which act like 
heat sinks (see Eq. (1/31)). Above 100 km, the heat input becomes zero, however 
the heat sink due to the advection flow remains. Therefore the temperature 
amplitude has a discontinuity there and its phase changes by 180°. This phe- 
nomenon is of course not to be ejected in reality since the heat input distribution 
is actually continuous. 

Vertical velocity and temperature have a constant time delay with respect 

to the solar heat input of Oh and 6h, respectively. The pressure changes 

sign at an altitude of z. =1.8 km in our model where z. .is the height 

isop isob 

of a true isobaric layer. Above that height, the pressure peaks at 18 h L. T. 
like the temperature, below that height it peaks at 6h L. T« 

The horizontal'velocity components are proportional to the pressure 
(see Eq. (4) and therefore have the same vertical distribution. 
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The latitudinal wind is shifted by 6h with respect to the pressure and peaks at 


Oh L. T. above n -' iso k* The longitudinal wind has the same phase as the pressure. 

A noon-midnight meridional cross-section of the global wind system of 
1 1/2 

the 0^’ -mode within the lower atmosphere has been drawn in Fig. 5. The 
wind rises at the zone of maximum heating on the day time hemisphere, crosses 
the poles and falls on the night time hemisphere. This wind system is com- 
pletely equivalent to the circulation cell of a thermally driven wind. The 
longitudinal winds (with maximum values at 18 h L. T. on the day time hemis- 
phere and at 6h L. T. on the night time hemisphere) are indicated in Fig. 7 
within the various sections of the cell by the symboles E = easterly (blowing 
from the east) and W = westerly (blowing from the west). 

In Fig. 6 we repeated these calculations for thermospheric conditions 
(using Eq (1/41)). Here we assumed a solar heat input of 


r 


for 

for 


z > z 


z < z 


(18) 


and adopted the temperature and dissipation data at 300 km altitude (see Fig. 2). 

The range z q < z < z q + 200 km in Fig. 6 represents the height region 

between about 200 and 400 km. As we already mentioned, the general behavior 
1 1/2 

of the 0j - mode at thermospheric heights is similar to the situation 
within the lower atmosphere as shown in Fig. 4. However, there exist some 
remarkable differences. First, the penetration depth of the quasi-evanescent 
mode is much greater than the corresponding value of Eq. (17) and of the order 
of 


Az ^1200 km . 
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Therefore, the assymptotic values of pressure and vertical wind indicated 

by the vertical dashed lines in Fig. 6 are still twice the values at z-z q = 200 km. 

/ 

Second, the pressure field has a minimum but final value near 4 km above the 
lower boundary of the heat input. Therefore, a true isobaric layer is not to 
be expected within the dissipative thermosphere. Moreover, the phase 
difference between the pressure below and above the quasi-isobaric layer is 
not greater than 10 h (Fig. 6b) as compared with exactly 12 h within the 
nondissipative lower atmosphere (Fig. 4). 

Within the first 200 km above z q the phases of pressure, temperature and 
vertical wind decrease slowly with altitude (Fig. 6b). Below the height z q 
which is the lower bour dary of the heat input only downgoing free internal 
waves exist which decay rapidly. This indicates the small influence the ther- 
mospheric region exerts on the lower atmosphere. 

The horizontal winds are proportional to the pressure and therefore have 
the same height distribution. The upper scale in Fig. 6a gives the magnitude of 
the lafitutdinal wind at the poles determined from Eqs. (4) and (10). The phase 
of the winds at the poles which can be derived from Figs, lb and lc are plotted 
in Fig. 6b and give phase delays with respect to the pressure of 11.5 h and 
5.5 h, respectively. Above 200 km the northerly wind peaks at about 2 h L. T. 
The westerly wund, shifted by 6 h, peaks at about 20 h L. T. 

The circulation cell due to the heat input above z q ~ 200 km is similar 
to the cell in F'g. 5. However, the change in sign of the longitudinal wind 
along one meridian, as indicated in Fig. 5 by the vertical lines at = 60°, dis- 
appears. Moreover, above the quasi-isobaric layer the maximum northerly 
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wind flows at 2 h L. T. as compared with Oh L.T. in Fig. 5. The total wind 

1 1/2 

circulation of the -mode is the sum of the individual circulation cells 

generated by the heat inputs within the various height regions including the 
regions between 100 and 200 km and above 400 km. If the heat input decreases 
nearly exponentially with height, we except only one fundamental circulation 
cell with the isobaric layer near the ground. Any individual cell above the 
fundamental cell can form only within a height range where the heat input has 
a relative maximum. This may happen within the Ozon layer or above the 
mesopause. 

The results of Figs. 6a and 6b justify our assumption about the height 
structure of the horizontal winds which was needed for the determination of the 
dissipation factors Z y j g in (1/3). 

3b Comparison with experimental data 

The numerical results shown in Figs. 1 to 6 are compared with 

available experimental data. From the work of Kato (1966), Stening (1969), 

Tarpley (1970) and Volland (1971) it appears to be well established that the 
1 1/2 

0^’ -mode is the main generator of the geomagnetic Sq current. Its 

horizontal wind field has the right latitudinal structure to drive the Sq current. 

1 l/2 

The phase of the wind of the 0^ ’ -mode changes slowly enough with altitude 

so that the induced electric current densities produce an optimally large 

magnetic field on the ground. From our result we may conclude furthermore 

that, provided the solar heat input peaks at local noon, the phase of the 

1 1/2 

latitudinal wind of the 0^’ -mode varies from 0 h L. T. below 100 km to 
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about 2 h L. T. above 300 km (see Figs. 4 and 6). This is just the phase value 

necessary in a consistent dynamo theory. Volland (1971) used a phase of 

u of 0 h L. T. which is the expected phase in the nondissipative atmosphere. 

Tarpley (1970) assumed a phase of u of 1. 3 h L. T. which is a value to be 

expected above 100 km. We moreover can estimate from the magnitude of the 

wind velocity necessary to explain the observed geomagnetic Sq variation (which 

is u ~ 120 m/sec) and from our result in Fig. 4 the normalized heat input of 
max 


J 


1 , 1/2 

1 




0.05 


(19) 


which is one order of magnitude greater than the corresponding values estimated 

for the Ozon-layer (Chapman and Lindzen, 1970). Finally it can be shown that 

the structure of the Sq current changes only slightly if the wind field of the 
1 1 /2 

0^ -mode changes with height according to Fig. lb and lc. 

We now compare our results of Figs. 1 to 6 with available satellite drag 
data. As it is obvious from Eqs. (1/9) and (1/13), the predominant diurnal 
component of density and temperature variations within the thermosphere 
correspond to the fundamental heat input component P*. Assuming a maximum 
heat input during local noon, the component of the density aiid temperature 
variation peaks with a time delay of about 3 h at 15 h L. T. Both latitudinal 
structure and time delay of the diurnal component are completely consistent 
with our result given in Figs. 1 and 6. 

The experimental facts , namely the tendency of an increasing relative 
pressure amplitude with height approaching an assymptotic value and the 
nearly isothermal behavior of the thermosphere are also indicated in Fig. 6. 
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Here, we should mention again that our temperature values appear to be 
underestimated. Especially, the slow decrease of the temperature amplitude 
with height is not realistic. . This results from the neglect of heat conduction 
waves which play an increasingly significant role at thermospheric heights 
(see the results of full wave calculation in paper II and the paper by Volland 
and Mayer (1970)). 

For a comparison between the relative pressure amplitude in Fig. 6 and 
corresponding abservational data (e.g., Jacchia 1965) we have to take into 
account that a pressure component generated by a heat input below z q ~ 200 km 
height must be added to the pressure amplitude in Fig. 6 to obtain the total 
observed pressure variation (Volland and Mayr, 1970). The total pressure 
amplitude at 400 km height may therefore have reached nearly the assymptotic 
value indicated by the vertical dashed line. Taking a total pressure amplitude 
at 400 km of about 2/3 of the assymtotic value in Fig. 6 and comparing it with 
the observed relative pressure amplitude of about 0. 5 at 400 km we find a 
normalized heat input of 

jJ’ 1/2 ~-0.1 (20) 

which is not inconsistent with the result of Eq. (19). 

The average height integrated heat input per area necessary to generate 
the observed density amplitude at the thermospheric heights can be estimated 
from Eqs. (1/9) and (1/27). It is for the diurnal component 


a 


2 Q 


4,1/2 


xuv 
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Thus it is 


Q XUV 

47Tr 2 
o 


= Q 


xuv 



— 'V 1 erg/cm 2 sec 
2 K 1 


-2 3 

above 120 km height with p ~2. 6 x 10 erg/cm at 120 km altitude, H ~40 km 

P ® 

1 1/2 

and J^’ from Eq. (20). This is a plausible number for the effective solar 

XUV heat input within a column above 120 km height (Hinteregger et al. , 1965). 

1 l/2 

Finally, it can be shown that the horizontal wind structure of the 9 ’ 

-mode at 300 km altitude is essentially identical with the wind fields derived 
from Jacchia’s pressure fields (Jacchia and Slowey 1967) by Geisler (1967) and 
Kohl and King (1967). 


4. THE ANTISYMMETRIC DIURNAL TIDAL MODE©*’ 1 / 2 
4a Eigenvalue and equivalent depth 

1 1/2 

The antisymmetric diurnal tidal mode 6 ’ which degenerates to the 
spherical function P* at thermospheric heights is according to Eqs. (1/9) 
and (1/13) the most important antisymmetric diurnal mode and describes 
the annual modulation of the diurnal variation of the solar heat input and 
its corresponding density variation which is proportional to 

0i’ l/2 cos T cosft t . 

2 a 

It will be shown in this section that within the lower nondissipative 
atmosphere that mode also degenerates to the spherical function P*. 

It is known from the work of Hough (1898) that the spherical function P * 


/ 
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must be added to the set of Hough functions to make the system complete. 

Therefore, P* indeed is itself a Hough function within the nondissipative 

atmosphere. In order to study its dependence on altitude, we proceed in a 

manner outlined in section 3a and start from the expressions 

= p 1 + 6 1 ’*- P 1 
2 2,4 4 

= - - . J — — — (ih 1 ’ f p 1 + pA 

*2 2f Cysind \ 2, 1 1 ^2,3 3 J 


(21) 


<t>: 


l,f 


4.f 




2 " 2f v ^ysint? \ 2* 2 P 2 +<P 2i2 P 4/ * 

This set of equations fulfilles again exactly the two first equations of 
(3) and leads to the relationship (with again f = f = f^) 

*2:! = - f k (4 - 4 /l2 - f k + 1/12 ) / < 5A) — - ^ 

*l'\ = 4 ^ f k ( f kin +1/4 ) ( f k 2 + V 2 - i/2 ) 7 < 5A > - 0 

^2,2 = f k (4 + 3f k /28 ' V 7 + 9/28 ) /A — - 1 
* 2 , 4 = 12 Vf f k( f k + V 2 ' 1/2 ) ' < 35A )— 0 

6 = - 3^10 (f£ + f k /2 - l/2j / (35A)-— 0 


(22) 


with 


A = f ^ + f^/4 - 27 f^/28 - 3 f k /28 + 3/28 — - 3/32 


where again the r. h. s. of (22) is determined for f k = f = 1/2. From (22) and 


(7) we find an equivalent depth of 


l,f _ 


2f0 2 r 2 A 


3g ( f k ' 1/12 ) + f k /2 ‘ 1/2 ) 


OO 


(23) 
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which becomes infinite for the diurnal component (f^ = f = l/ 2 ) within the 


nondissipative atmosphere. 

As it is clear from Eqs. (3), (7) and ( 22 ), our solution for the Hough 


function is exact within the nondissipative atmosphere and gives 


11/2 1 

0 9 ’ — ► P . Within the dissipative atmosphere it is for 


f. 


> 1 


6 1,f «l/ff 
2.2 ‘ k 


and 


<P 


l,f 

2,4 


l/l 


A k 




0 


l,f 

2,2 

l.f 

2,1 


-3 V3/5 : 


4 V^/5 


(24) 


approaching the formulae (1/20) and (1/26) (Note that 


dP 


d# 


{ * p 5 • » p !) 


4b Wind and pressure fields 

In order to study the behavior oi the -mode within the atmosphere 

we did the same calculations as in section 3b and plotted in Fig. 7a the 
1 1/2 

eigenfunction O ’ versus co-latitude for four different heights in amplitude 

and phase. Here we note a significant deviation from the structure of P 9 only 

within the height range between 100 and 200 km. Figs. 7b and 7c give the 
ll/2 ll/2 

wind function and ^9 versus co-latitude for the four heights 

considered. Here too a transition region exists between 1 00 and 200 km. Above 
and below this region the structure of the longitudinal wind is the same and 
proportional to 

Pg / sin# = \/3 cos# 
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However the structure of the latitudinal wind changes from 

P* / sini? = 1 

within the lower atmosphere to a structure proportional to 

dP* = VT (2 cos 2 d - l) 
d# 

within the upper part of the thermosphere. 

The propogation factor O' and the attenuation factor (3 - X of the ©;* * 

-mode are plotted in Fig. 3. We observe a strong increase of the attenuation 

factor between 100 and 200 km and then a slow decrease with height. 

The attenuation factor $ - 1 is greater by a factor of about 3 than the 

i 1/2 

corresponding value of the 0^’ -mode. 

1 1/2 

We also did calculations of the generation and propogation of the O' 

-mode within the lower atmpsphere (Fig. 8) and within the thermosphere 
(Fig. 9) adopting the same assumptions as in section 3b. We assumed a solar 
heat input of 

jl’ 1/2 = - 0.01 

below 100 km in Fig. 8 and above z q ~200 km in Fig 9 which is equivalent 
to maximum heating at local noon during northern summer (June -solstice). 

In the case of the nondissipative atmosphere our solution of the eigenvalue q 
(1/37) degenerates to 

q l. i/ 2 = - jA (25) 
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1 1 / 2 

because of h 2 ’ ' —►<». Thus, we find from (1/50) 

F = - 2jA 
a J 

F b - 2jA(y - 1) 

r =- J 
a A(e - 1) 

G, = 0 
b 


( 26 ) 


which means , no downgoing wave is generated. For the upgoing wave we find 
with € = 1: 


and 


lim a(z) = k J z 
C-l ° 

w/c o = k Q J Z 
p/p o = -2jAk Q J z 
T/T q = -jJ 


(27) 


Within the lower atmosphere , this mode behaves like a degenerated 

evanescent mode with constant phase . however with the relative amplitude 

of the pressure linearily increasing with height (see Fig. 8). 

11/2 

Within the termosphere (Fig. 9) the © 2 ’ -mode behaves like a quasi- 

1 1/2 

evanescent mode and is similar to the 0^* -mode in Fig. 6. However, 

there exists one striking difference. As we note from Fig. 3, the attenuation 

1 1 /2 

factor # - 1 of the 0 9 ’ -mode is much greater than the corresponding value 
1 1 /2 

of the 0j “ -mode Therefore the penetration depth Az (see Equ. (17)) is 
much smaller, and the assymptotic values of vertical wind and pressure 
(indicated by vertical bars in Fig. 9a) are reached at a lower height than 
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in the case of Fig. 6a. On the other hand, the magnitudes of the assymptotic 
values are not affected by the attenuation factors (see Eq. (1/47) and the 
vertical dashed lines). 

The circulation cells corresponding to the heat inputs within the lower 
atmosphere (Fig. 8) and within the thermosphere (Fig. 9) are shown in 
Figs. 10a and 10b. Both ceils differ significantly from each other. Within 
the lower atmosphere the isobaric layer is the ground. Thus , only one cell 
extended over the entire globe is formed, no return wind exists and the lati- 
tudinal wind blows from south to north during northern summer and noon. A 
westerly wind (directed toward the east) blows on the northern day time hemi- 
sphere at 18 h L. T. Within the thermosphere (Fig. 10b), two cells on each 
hemisphere are formed separated from each other at latitudes ±45°. These 
winds change their direction during northern winter. 


4c Comparison with experimental data 

We now want to compare the results of the forgoing sections with available 

data. From an analysis of the geomagnetic Sq variations, Volland (1971) 

derived an antisymmetric wind field at E -layer heights which has the exact 

1 1/2 

form of the wind of the tidal 6 ’ -mode given in Figs. 7b and 7c at z = 100 
km. The wind velocity at the poles was determined to be 


u 


1 , 1/2 

2 


20 m/sec 


The phase of the latitudinal wind during northern summer was 0 h L. T. 

. 1 , 1/2 


which agrees with the phase of u f 


plotted in Figs. 8 and 10a if the heat 
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input peaks at local noon. The ratio between the antisymmetric and the 
symmetric diurnal winds derived from the Sq current was 


u 


u 


1 , 1/2 

2 

1 , 1/2 


pi' 1/2 

1 , 1/2 

P 1 


20 

V3 105 


= 0.11 


d=0 


( 28 ) 


From the Jacchia model (Jacchia and Slowey 1964) we find a (2,1)- 
component of the exospheric temperature (Eq. (1/14)) of 


T 2 P 2 0. 008 



0. 07 


(29) 


and a phase value for T 2 (that is: of the pressure component p*) of 15 h L. T.* 
This phase agrees sufficiently well with the value of 14 h L. T. in our 
theory (see Fig. 9b). We compare the amplitude ratio of (29) with the 
corresponding value in our theory which gives 


t 1,1/2 1,1/2 -1,1/2 

13 = £2 _ 1 , j 2 , . _ J_ .5. V T. 5 _ o 14 

1072 1,1/2 3 -1,1/2 3 2x8 0 

T 1 Pi J 1 


(30) 


*Note that Jacchia’s exospheric temperature represents merely the density 
amplitude above about 200 km. A relative temperature amplitude in the 
Jacchia-mode of T^/T q = 0. 1 corresponds to a relative density amplitude 
of pJJ 3 / P Q ^ 0.4. Moreover, the density amplitude is proportional to the 
pressure amplitude, and it is P^ 1 / P Q ~ 0. 85 P^/p o within thermospheric 
heights above 200 km (see also Equ. (1/45)). 
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Here, for the heat input components we use the XUV-values in (1/9), 

(1/27) and (1/47). This number is greater by a factor of two than the experi- 
mental value in (29) but it appears to be sonsistent with the number of (28) 

Considering this discrepancy, we should take into account tha f the 
Jacchia-model is based on a plausible but arbitrary temperature and there- 
fore pressure profile which is fixed for all possible frequencies of the 
disturbance, which is our notation means, fixed for all wave domain numbers 
(n, m, f). As we may notice from Figs. 6a and 9a, this is certainly not the 
case. The temperature and pressure profiles strongly depend on the wave 
domain numbers. Moreover, the exospheric temperature distribution of the 
Jacchia-model is described by an arbitrary analytic function which appears 
to fit well the main c ^mponents like the symmetric diurnal (1, l)-component, 
but which may more or less deviate from second order components like 
the antisymmetric diurnal (2, 1) -component. 

On the other hand, the heat input distribution in Eq. (1/9) based on rather 
crude assumptions as well as our simplified theory can of course give only 
first order approximations. In view of these uncertainties, the numbers given 
in (28), (29) and (30) appear to be quite consistent with each other agreeing 
completely in phase and leading to an amplitude ratio between the anti- 
symmetric and the symmetric diurnal component of the order of 0. 1. 
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5. SYMMETRIC SEMIDIURNAL TIDAL MODE0 l’ 1 

£ 


5a Eigenfunction and equivalent depth 

2 1 

The symmetric semidiurnal tidal mode 0 ’ degenerates within the 

& 

o 

thermosphere to the spherical function P„ and behaves like a quasi -evanescent 
mode there. It is responsible for the observed semidiurnal component of the 
density in Eq. (1/14). Within the lower atmosphere, this mode is a propaga- 
tion mode called the (2,2) tidal mode in the literature (e. g. , Chapman and 
Lindzen, 1970). It has been tabulated by Chapman and Lindzen and has the 
following form: 

0(2,2) = 1.08 Pg - 0.368 P^ + 0. 044 Pg + . . . (31) 


where the neglected terms are smaller than 3% of the first term. In order to 
study its height dependence within the nondissipative atmosphere we start 
in the usual way with the expressions 


0 2 ’ f = p 2 + 6 2 ’^ P 2 
2 2 2,4 4 


* 2 > f = J -f’f 

2 2f STsind ^2,3 
u 



<t> 


2,f 

2 


2 

2f v £ysind 





(32) 


and obtain from the first two equations in (3) the exact relationship (again 


wlth f u = f v f k> 
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= 2 f k (f k + 1/2)/ (vs^)— - - 7/vr 

*2 . 2 =f k( f k- 5/14 ) /A — 3/2 
<^4 = 2vU f k /(35A) ► 2/VI5 

6 2 ’^ = - Vl 5 /( 35 ^ ► - 1/ Vl5 

with 

A = - f k /2 - 1/14 — — 3/7 


From (32) and (7) we obtain the equivalent depth of 


,2, f 2fH 2 r 2 A 

2 3g (f k - 1/6) 


12 Q 2 r 2 
35 g 


7.47 km 


(33) 


(34) 


In (32) and (33), the r. h. s. numbers are valid for f^ = f = 1. 

We can compare the exact numerical value of the (2,2)-mode which is 

7. 85 km (Chapman and Lindzen, 1970) with our approximate value in (33) and 

find an agreement within about 5 %. The bottom curves in Fig. 11a show the 

exact latitudinal structure of the (2.2)-mode (dash-dotted line) and our 

approximated function (31) with the parameters of the r. h. s. of (32) (solid line) 

and again reveal satisfactory agreement. The bottom curves in Figs, lib and 

2 1 2 1 

11c give the latitudinal and the longitudinal winds functions and ^2 

according to our approximation (solid lines) and according to exact theory for 
2 1 

the © 2 * -mode (dash-dotted lines) within the nondissipative atmosphere. Here 
too the agreement is sufficient, although not as good as for the eigenfunction 
© 2 ’ 1 itself. 


28 


Within the dissipative atmosphere it is for fjl > 1 


>2 * f 

2,4 


,„2,f 

(n 7 oz 

V 2,4 


1/it 


VI *2.f 

2 *2,3 


<P 


2, f 

2,2 


1/f 


K 


1 , 


and the wave amplitudes of pressure and winds again approach the forms of 
Eqs. (1/20) and (1/26) degenerating to the ordinary spherical functions. 
(Note that 


d# y/K sin# 


5b. Wind and pressure fields 

After repeating the calculations done in sections 3b and 4b we plotted in 

2 1 

Fig. 11a the eigenfunction © 2 ’ versus co -latitude for the four different 

2 1 

heights considered. Above 200 km and below 100 km, the © 2 ’ -mode has 
nearly the same structure with constant phase independent of #. At 150 km, 
however, we notice a significant phase transition with co-latitude. 

In Figs lib and 11c the wind function $ 2 ’ an< * 4> 2 ’ are plotted versus 
co-latitude for the four different heights. Here the structure of the latitudinal 
wind (Fig. lib) does not change with height in our approximation and remains 
proportional to 


3 

sin0 


■v? 


sin# cos# 
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However, the structure of the longitudinal wind changes significantly with 
height approaching the ^-dependence of 



The propagation factor ck and the attenuation factor $-1 of the 0^’ 1 -mode 

are plotted in Fig 3. Here, the propagation factor is different from zero at 

and below 100 km indicating a behavior for this mode resembling a propagation 

mode. Likewise, the attenuation factor - 1 is negative below 165 km (the 

dashed part of the curve in Fig. 3a), corresponding to a net exponential 

increase with height of the relative wave amplitudes. However, above 165 km 
2 1 

the © 2 ’ -mode has changed into a quasi-evanescent mode and behaves 
essentially like the other modes considered. 

Simulating lower atmospheric (Fig. 12) and thermospheric (Fig. 13) 
conditions we calculated the wave amplitudes versus heights due to a heat 
input of 

9 1 

Jg’ = 0.01 

which peaks at local noon and at local midnight. Since within the lower 

atmosphere of scale height H q = 6 km (the value used in Figs. 4 and 8), 

2 1 

the 0 ’ -mode is already a trapped mode, and since we want for illustration 
to discuss a real propagation mode, we applied a pressure scale height of 
H q = 7. 8 km in Fig. 12 corresponding to a mean temperature of T q = 260° 

The result of Fig. 12a clearly reveals the characteristics of a propagation 
mode, namely the exponential increase of the relative wave amplitudes with 
height apart from the pressure amplitude near the ground where reflection 
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properties have to be taken into account. Moreover, the phase in Fig. 12b 
decreases with altitude as we expect from an upward propagating gravity wave. 
The free internal propagation above 100 km— the termination of the heat input 
in our model— does not alter that propagation behavior. 

The relative amplitudes in Fig. 12a can not increase indefinitely with 

height in reality. If the relative wave amplitudes approach the value one, 

our perturbation theory breaks down, and the wave mode will transfer into 

a shock wave with wave energy dissipation at the shock front. Moreover, 

higher harmonics of the semidiurnal wave may be generated. On the other 

hand, within the critical height above 100 km where such transfer may happen, 

2 1 

the propagation properties of the 0 ’ -mode change drastically (see Fig. 3) 
which prevents furthermore an unlimited increase of the wave amplitude. 

This transition from a propagation mode into a quasi-evanescent mode must 
be treated by numerical full wave calculations and will be done in paper II. 
Within the thermosphere (Fig. 13), our wave mode behaves quite normally 
like a quasi-evanescent mode. 

2 1 

The meridional cross-section of the wind system of the © ’ -mode due 
to the assumed heat input within the thermosphere is given in Fig. 14. Two 
circulation cells are built up. Maximum latitudinal winds directed to the 
poles blow at about 7. 5 and 19.5 h L. T. above the quasi-isobaric layer ^ 
The longitudinal winds, shifted by 2. 5 to 3 hours with respect to the latitudinal 
winds, are westerly at 5 and 17 h L. T. at the equator and above z. ^ and have 
the same direction over the entire day time hemisphere. 
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The V7 s nd structure within the lower atmosphere is similar to that of 
Fig. 14 except for the phase. Here (Fig. 12b), the phase changes continuously 
with height reaching a phase transition of 180° after about 100 km. Contrary to 
that behavior , within the thermosphere such transition occurs within a height 
range of about 20 km (Fig. 13b). 

5c Comparison with experimental data 

In comparing the result of the forgoing sections with experimental data, 

we have to take into account that the primary source of the semi-diurnal tidal 

component is absorption of Ozon in the stratosphere and mesosphere (Bulter 

and Small, 1963) which means a more complicated vertical structure of the heat 

input than assumed in Fig. 12. Second, due to the nearly equal values of 

2 1 

the pressure scale height and of the equivalent depth of the © 2 ’ -mode within 

the lower atmosphere , the real mean temperature profile is important for the 

generation and propagation of that mode. E. g. , within the mesosphere this 

mode is evanescent. Therefore, the results in Fig. 12 which were based on an 

2 1 

isothermal atmosphere and a constant relative heat input J ’ served merely for 

illustration and cannot be used for a quantitative comparison. More realistic 

model calculations of the fundamental semidiurnal mode may be found in the 

book of Chapman and Lindzen (1970). 

2 1 

Since the © 2 ’ -mode is partly trapped within the lower atmosphere, the 

2 1 

(2,4)-mode (in our notation: the © 4 ’ -mode) which is a pure propagation 
mode throughout the lower atmpsphere may significantly contribute at 
thermospheric heights. The wind data obtained between 80 and 100 km from 
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Meteor trail measurements (Elford, 1959 ; Greenhow and Neufield, 1961) which 

give velocities of the order of 20 m/sec for the semidiurnal component there - 

2 1 2 1 

fore may be related to theG^’ -mode rather than to the 0 9 ’ -mode (Hines, 

1963). Similarly, the evaluation of a sequence of wind measurements at 

intervals throughout a single night above about 120 km (Murphy and Bull, 1968) 

leading to vertical wave lengths of 50 - 60 km give further support for the 

2, 1 

predominance of the G^ -mode at those heights . 

On the other hand, within the upper part of the thermosphere the© 2, 1 

4 

-mode should be suppressed due to its quasi-evanescent behavior like (see 
Eq. 1/47) 


2,1 


,2, 1 

P 4 

^ 1 

J 4 

1,1/2 

P 1 

~ 10 

,1,1/2 

J 1 


1 

10 


(35) 


and should become insignificant there. This idea is consistent with the 

density da\;a deduced from satellite drag measurements (see Wq. (1/14)). 

2 1 

However, the semidiurnal G^’ -mode should be related to the fundamental 
diurnal ’ -mode by 


2,1 

P 2 

^ 1 

t 2,1 

J 2 

1 

5 \/l2 

1,1/2 

P 1 

~ 3 

,1,1/2 

J 1 

~ 3 

2 x 16 


(36) 


if we adopt the numbers given in Eqs. (1/9), (1/27) and (1/47), and the time of 

2 1 

maximum of the pressure component of the0 2 ’ -mode is about 13 - 13. 5 h. 

L. T. according to Fig. 13b. 
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We compare these numbers with the experimental result of Eq. (i/l) 


which gives 



0. 035 
0.12 


0.29 


(37) 


and a phase of 13.2 h. L. T. . 

While the observed phase value excellently agrees with the theory, the 

observed amplitude ratio in (37) is greater by about 50% than the calculated 

number in (36), This difference can be interpreted as resulting from a 

relatively strong influence of the wave generated below 200 km where the 
2 1 

0g’ -mode transfers from a propagation mode into a quasi-evanescent 
mode. Indeed, Lindzen (1970) calculated that the semidiurnal component 
from below should be significant at thermospheric heights although according 
to our theory his result appears to be over-estimated. A final answer to that 
problem can be given only in a numerical full wave treatment (to be done in 
paper II) where the transition range between 100 and 200 km must properly be 
taken into account. Here also the question may be solved why the phase of the 
semidiurnal component of the neutral temperature derived from Thompson 
backscatter measurements at 250 km (Mahajan, 1969; Waldteufel and McClure, 
1969) is phase shifted by about 2 hours with respect to the corresponding 
maximum of the density (Volland, 1970). 


6. CONCLUSION 

Based upon a theoretical treatment of a threedimensional spherical 
thermospheric model outlined in a previous paper (Volland and Mayr, 1971a) 
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we determined the eigenfunctions and the equivalent depths of three fundamental 
tidal modes— the symmetric diurnal mode, the antisymmetric diurnal mode and 
the symmetric semidiurnal mode— and calculated their dependence on height 
and latitude within the dissipative thermosphere. The corresponding pressure 
and horizontal wind fields of these modes have been calculated and meridional 
cross-sections of the circulation cells due to these modes have been determined. 
The theoretical magnitude and phase of the wave amplitudes of pressure and 
horizontal winds of the various tidal wave modes have been compared with 
corresponding available experimental data and show satisfying agreement. We 
conclude from this results that a spatial solar XUV heat input distribution 
within the thermosphere proportional to the solar zenith angle is the main 
generator of the three tidal waves considered and that the Jacchia -model of 
the thermospheric density distribution (Jacchia and Slowey, 1967) can readily 
be interpreted as the thermospheric density response to that energy input. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Dependence on co-latitude t?of the eigenfunction of the symmetric 
1 l/2 

diurnal tidal 0^ * -mode for four different heights 

1 l/2 

Fig. la: Pressure function 0^’ 

1 l/2 

Fig. lb: Latitudinal wind function 

1 l/2 

Fig. 1c: Longitudinal wind function 

Solid lines: magnitudes 

Dashed lines: time delay (in hours) 

Dash-dotted lines: exact values within the nondissipative 
atmosphere 

Vertical profiles of mean temperature T q and dissipation factors 
due to heat conduction (Z^), and due to viscosity and ion drag 
(Z^in) used in our model calculations. 

Fig. 3a: Attenuation factors /3 —1 ( j3 is the imaginary part of 

2 1 

the eigenvalue). $ ^ -1 is negative below 165 km (the dashed 

line in Fig. 3a), 

Fig. 3b: Propagation factors a ( a is the real part of the 
eigenvalue) . 

Relative wave amplitudes of pressure p, temperature T and 

1 1/2 

vertical wind w versus altitude of the symmetric diurnal 0^* 
-mode within the lower nondissipative atmosphere. The latitudinal 
velocity u is proportional to the pressure. Its pole values are 
given in the upper right scale. The time delay r Q gives the time 
of maximum amplitude in local time. The heat peaks at local 
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noon and is terminated at 100 km. Therefore above 100 km only 
upward propagating free internal wave exist. The pressure 


Figure 5. 


Figure 6. 


Figure 7. 


amplitude is zero at the height of the isobaric layer z^ g ^ =1.8 km 

and shifts in phase by 12 hours there (see the box on the l.h.s.) 

Noon -midnight meridional cross-section of the wind system of the 
1 l/2 

tidal 0^ -mode within the lower atmosphere. The 
longitudinal winds are phase shifted by 6 hours with respect to 
the latitudinal winds. Maximum westerly winds (winds blowing 
from the west) are indicated by the symbol n W", easterly winds 
(winds blowing from the east) are indicated by ,, E M . 

Relative wave amplitudes of pressure p, temperature T and 


vertical wind w versus altitude of the tidal 0 


1, 1/2 .. 


mode within 


the dissipative thermosphere 200 kmj. The generating heat 

input peaks at local noon and is terminated at the height z q . Below 

that height only downgoing free internal waves exist. 

Fig. 6a: magnitudes (the vertical dashed lines give the maximum 

asymptotic values of pressure and vertical wind). 

Fig. 6b: times of maximum (in hours local time). 

Dependence on co -latitude d of the eigenfunction of the 

1 l/2 

antisymmetric diurnal tidal 0 2 * -mode for four different 


heights. 


Fig. 7a: Pressure function 0 2 


1, 1/2 


1, 1/2 


Fig. 7b: Latitudinal wind function 

1 l/2 

Fig. 7c: Longitudinal wind function 4> 2 * 
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Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


Solid lines: magnitude 

Dashed lines: time delay (in hours) 

Relative wave amplitudes of pressure p, temperature T and 

vertical wind w versus altitude of the antisymmetric diurnal tidal 

® 2* ^ ” m0( * e within the lower nondissipative atmosphere and 

during June-solstice. The generating heat input peaks at local 

noon and is terminated at 100 km. The time delay r gives the 

time of maximum amplitude in local time. 

Relative wave amplitudes of pressure p, temperature T and 

1 l/2 

vertical wind w versus altitude of the tidal 0 2 * -mode within 
the dissipative thermosphere during June-solstice. The generating 
heat input peaks at local noon and is terminated at z q -w 200 km. 

Fig. 9a: magnitudes (the vertical dashed lines give the maximum 
assymptotic values of pressure and vertical wind). 

Fig. 9b: times of maximum (in hours local time). 

1 l/2 

Meridional cross-section of the wind system of the tidal 0 * 

Li 

-mode during June-solstice within the lower atmosphere and within 

the thermosphere. Maximum longitudinal winds phase shifted by 

about 6 h with respect to the meridional winds are indicated by 

"W" (westerly winds) and by M E M (easterly winds). 

Dependence on co-latitude d of the eigenfunction of the symmetric 

2 1 

semidiurnal tidal 0 2 ’ -mode for four different heights. 

2 1 

Fig. 11a: Pressure function 0 2 * 

2 1 

Fig. lib: Latitudinal wind function ♦ 2 * 
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2 1 

Fig. 11c: Longitudinal wind function 4> 2 * 

Solid lines: magnitudes 

Dashed lines: time delay (in hours) 

Dash-dotted lines: exact values within the nondissipative lower 
atmosphere. 

Figure 12. Relative wave amplitude in magnitude (Fig. 12a) and time of 

maximum in hours L. T. (Fig. 12b) of pressure p, temperature 

T and vertical wind w versus altitude of the symmetric 

2 1 

semidiurnal tidal 0 2 * -mode within the lower nondissipative 
atmosphere. The longitudinal wind v is proportional to the 
pressure. Its equator values are given in the upper right scale. 

The generating solar heat input peaks at local noon and midnight 
and is terminated at 100 km. 

Figure 13. Relative wave amplitudes in magnitude (Fig. 13a) and time of 

maximum in hours local time (Fig. 13b) of pressure p, temperature 

2 1 

T and vertical wind w versus altitude of the tidal 0 2 * -mode 

within the dissipative thermosphere. The generating solar heat 

input peaks at local noon and midnight and is terminated at 

z ~ 200 km. 
o 

2 1 

Figure 14. Meridional cross-section of the wind system of the tidal 0 2 * 

-mode within the thermosphere. Maximum longitudinal winds phase 
shifted by about 3 hours with respect to the latitudinal winds are 
indicated by ,f W ,f (westerly winds) and by n E ,f (easterly winds). 


42 











ATTENUATION 



Fig. 3d Attenuation factors )3-l (0 is the imaginary part of the eigenvalue). /J?* 1 -1 is negative below 
165 km (the dashed line in Fig. 3a). 

Fig. 3b Propagation factors a (a is the reel part of the eigenvalue). 
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that height only downgoing free internal waves exist. 

Fig. 6a: magnitudes (the vertical dashed lines give the maximum asymptotic values of pressure and vertical wind). 
Fig. 6b: times of maximum (in hours local time). 
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Figure 7. Dependence on co-latitude i? of the eigenfunction of the antisymmetric diurnal tidal @ 2 '* 2 -mode for four different heights. 
Fig. 7a: Pressure function © 2 ' “ Fig. 7b: Latitudinal wind function Fig. 7c: Longitudinal wind 

. . function 3>i ,1/2 

Solid linfco, magnitude Dashed lines: time delay (in hours) 
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WAVES AMPLITUDES OF 0$ 

, 10 4 * w/c„ 


MODE 


p 

(18 h LT) 


P/Po and T/T 0 


Relative wave amplitudes of pressure p, temperature T and vertical wind w versus altitude 
of the antisymmetric diurnal tidal @ 2 ,1/2 -mode within the lower riondissipativc -itmos- 
phere and during June-solstice. The generating heat input peaks at local noon and is 
terminated at 100 km. The time delay r Q gives the time of maximum amplitude in local 
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Figure 12. Relative wave amplitude in magnitude (Fig. 12a) and time of maximum in hours L. T (Fig. 12h l i of pressure p, temperature T 
and vertical wind w versus altitude of the symmetric semidiurnal tidal mode \v ; i the lower nondissipative atmosphere. 
The longitudinal wind v is proportional to the pressure. Its equator values are given in the upper right scale. The gemmating 
solar heat input peaks at local noon and midnight and is terminated at 100 km. . 
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